The phase structure of QGP-Hadron in a statistical model using Cornell, Richardson 

and Peshier potentials 
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We study the phase structure of the QGP-Hadron system under quasi-static equilibrium using the 
Ramanathan et al. statistical model for the QGP fireball formation in a hadronic medium. While in 
the earlier published studies we had used the Peshier effective potential which is appropriate for the 
deconfined QGP phase but could be extrapolated to the transition region from the higher momentum 
regime, in this paper we study the same system using the Cornell and Richardson potentials which 
are more relevant for the low momentum confinement regime, but could again be extrapolated to 
the transition region from below. Surprisingly, the overall picture in both the cases are quite similar 
with minor divergences, (though,the results with the Richardson potential shows a sizable deviation 
from the other two potentials), thus indicating the robustness of the model and its selfconsistency. 
The result of our numerical results pertaining to the variation of the velocity of sound in the 
QGP-Hadron medium with temperature in the various scenarios considered by us, is that, the 
phase transition seems to be a gentle roll-over of phases rather than a sharp transition of either the 
first or second order, a result in conformity with recent lattice calculations, but with much less effort. 
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I. INTRODUCTION 

The formation of QGP droplet (fireball) is one 
of the most exciting possibility in ultra relativis- 
tic heavy ion collisions(URHlC) [l]The physics 
of such an event is very complicated and to ex- 
tract meaningful results from a rigorous use of 
QCD appropriate to this physical system is al- 
most intractable though heroic efforts at lat- 
tice estimation of the problem has been going 
on for quite some time [2]. One way out is 
to replicate the approximation schemes which 
have served as theoretical tools in understanding 
equally complicated atomic and nuclear systems 
in atomic and nuclear physics in the context of 
QGP droplet formation. This approach lays no 
claim to rigour or ab-initio understanding of the 
phenomenon but lays the framework on which 
more rigorous structures may be built depending 
on the phenomenological success of the model as 
and when testable data emerges from ongoing 
experiments. 

The nucleation process is driven by statistical 
fluctuations which produce the QGP droplets in 
a hadronic phase,the size of the fluctuations be- 



ing determined by the critical free energy differ- 
ence between the two phases. The Kapusta et al 
model (3) uses the liquid drop model expansion 
for this, as given by 



AF = —R^[P^^,(T,fiB) - P,,giT,fiB)] 



+A'KR^a + TcritT In 



q,g 



The first term represents the volume contribu- 
tion,the second term is the surface contribution 
where a is the surface tension, and the last term 
is the so called shape contribution. The shape 
contribution is an entropy term on account of 
fluctuations in droplet shape which we may ig- 
nore in the lowest order approximation. The crit- 
ical radius Rc can be obtained by minimising (1) 
withrespect to the droplet radius R,which in the 
Linde approximation [6] is. 
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II. THE STATISTICAL MODEL 
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In the approximation scheme of the Ra- 
manathan et.al [4 ] the relativistic density of 
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states for the quarks and gluons is constructed 
adapting the procedures of the Thomas-Fermi 
construction of the electronic density of states 
for complex atoms and the Bethe density of 
states[6] for nucleons in complex nuclei as tem- 
plates .The expression for the density of states 
for the quarks and gluons(q,g) in this model is: 



With the surface free-energy given by a modi- 
fied Weyl[7] expression: 

Finterface = '^^ j dkpy,eyl{k)5{k - T) , (8) 

where the hydrodynamical flow parameter for 
the surface is: 



(3) 

where k is the relativistic four-momentum of 
the quarks and the ghions, v is the volume of 
the fire ball taken to be a constant in the first 
approximation and V is a suitable confining po- 
tential relevant to the current quarks and gluons 
in the QGP [4] given by: 

V^^{k)k' = (1/2A;)7,,, g\k)T' - ml/2k . (4) 

where m is the mass of the quark which we 
take as zero for the up and down quarks and 
150 MeV for the strange quarks.The g(k) is the 
QCD running coupling constant given by 

g\k) = (4/3)(127r/27){l/ln(l + A;7A^)} , (5) 

where A is the QCD scale taken to be 
150 MeV In fact, eq.(3) may be considered as 
the starting ansatz of the model without further 
elaboration as far as our calculational program 
is concerned. 

The model has a natural low energy cut off at: 

(6) 

with N = [(4/3)(127r/27)]3. 

The free energy of the respective cases 
i(quarks,gluons,interface etc.) for Fermions and 
Bosons(upper sign or lower sign) can be com- 
puted using the following expression: 

F.^TTgiJ dkp,{k)\n{l±e-^V^y^) , (7) 



7 = X y/(l7y^Ta77^, (9) 

For the pion which for simplicity represents the 
hadronic medium in which the fire ball resides 
,the free energy is: 

= (3r/27r')i/ r^dklnil - g-V^^/^) . 
Jo 

(10) 

With these ingredients we can compute the 
free-energy change with respect to both the 
droplet radius and temperature to get a physical 
picture of the fireball formation,the nucleation 
rate governing the droplet formation, the nature 
of the phase transition etc. This can be done 
over a whole range of flow-paramctcr values [4], 
We exhibit only the two most promising scenar- 
ios in fig.l and fig. 2. Further investigating the 
properties of the corresponding free energies as 
in [4], it was found that the QGP-hadron phase 
transition is a weakly first order one. The figs, 
indicate that the set of flow parameter values 
used in fig. 1 lead to more stable droplets for 
a whole range of temperatures with experimen- 
tally significant fireballs of radius of the order 
of 5 fermi. The increase in the barrier height 
with temperature is indicative of the fact that 
the droplets will be more stable at higher tem- 
peratures. However, as the rate of droplet for- 
mation is controlled by the ratio of the barrier 
height to temperature, the probability of hav- 
ing more stable droplets is also smaller at higher 
temperatures. The S Vs. T graphs [4] by their 
smoothness indicate that the transition for a fi- 
nite chemical potential is neither a first order nor 
a second order one, but a roll over of phase from 
the hadronic to the QGP. 

We could have used any potential form in our 
model instead of the modified Peshier potential 
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(4) used above, the main reason for our use of 
this form was that the dynamical nature of the 
quarks in the QGP is well brought out by this. 
However, there is also a widespred use of the Cor- 
nell and Richardson potentials in the recent liter- 
ature [9, 10] in the context of QGP, though there 
is also misgiving in some quarters about the va- 
lidity of using these potentials defined mainly 
for static constituent quarks in the QGP sys- 
tem. For the sake of testing the robustness of our 
statistical model we have done the calculations 
using the Cornell potential [9] and Richardson 
Potential [9a] in our formalism. 

For Gluon Sector, the Cornell potential in mo- 
mentum space is 



where as=Strong coupling constant =0.09, 
C = Casimir Scaling = 9/4, a = String tension 
= 8T2, Tc = 170 MeV For Quark Sector, the 
Cornell potential for three flavors in momentum 
space is 



Vik) 



IGtto;, Stto" 
+ 



3P 



(12) 



where a,= 0.009, a = IST^, = 170 MeV 
For the Richardson potential in the momen- 
tum frame we use the form [9a] 

V{k) = (127r/27)(4/3)(lA2)(l/ln(l + k^A^)) 

(13) 

where A = 150 MeV and k is four momentum. 

We can thus redo our calculations using the 
above potential forms in our formulation, using 
the parameter values widely used in the htera- 
ture. 

III. RESULTS AND CONCLUSION 

The numerical results of our calculations is 
summarized by the graphs figs.l to 6, where the 
temperature variation of the velocity of sound 
in the QGP-Hadron system is displayed for the 
typical parameter values of the model for the 



200 220 240 200 

T (Mbv) 



FIG. 1: Variation of with temperature T at 
1g = Glq , Iq = 1/6 and /i = 300 with Peshier 
Potential . 
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FIG. 2: Variation of with temperature T with 
Cornell Potential. 



three different potentials namely the Peshier, the 
Cornell and the Richardson potentials. 

When we compare the numerical results aris- 
ing from the two different potentials we have con- 
sidered, the nature of phase transition remains 
the same in both the potentials' scenarios. The 
velocity of sound, whose square is the ratio of the 
entropy and specific heat at constant volume, for 
the two potentials are not too different as indi- 
cated by the graphs (1, 2). While the peshier 
potential gives almost a contant velocity over a 
whole range of temperature variation, so does 
the Cornell potential, except for a distinct kink 
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FIG. 3: Variation of with temperature T at 
7g = 6jq , 7g = 1/6 and A* = 0, 300 with Cornell 
Potential for only three flavored quark sector. 
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FIG. 4: Variation of with temperature T at 
7s ~ ^Oq ) 7q = 1/6 and = 0, 300 with Cornell 
Potential for pure gauge field. 

at around 160 MeV, the expected critical tem- 
perature for the transition. However, in the case 
of the Richardson Potential, the fall in Cg Vs. 
T is steeper. So, this difference in their predic- 
tion could be used to discard one of the potential 
forms which disagrees with the actual measure- 
ment of sound velocity in a QGP fireball for- 
mation event in the foreseeble future.The shift 
in the absolute values of the velocity in the two 
cases should not hold much physical significance, 
as the shift seems to be the result of some overall 
factor in the potentials used. As we have already 
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FIG. 5: Variation of with temperature T with 
Richardson potetial at jg = 6jq. 
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FIG. 6: Variation of Cg with temperature T with 
Richardson Potential at jg — 8jq. 

noted that we could have used any effective po- 
tential available in the literature of quark-gluon 
interaction, our present excercise only strength- 
ens our confidence in the robustness of the sta- 
tistical model used as two totally unrelated ef- 
fective potentials, one, the peshier, designed for 
the deconfined QGP phase and the Cornell for 
the confined lower momentum regime lead us ti 
similar scenarios in the QGP-Hadron transition 
region. Evidently the ideal gas hmit of C| = | 
is obeyed only by the graph in fig.l, indicating 
that the Cornell and Richardson potentials lead 
to violation of this limit, adding to our already 
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stated reasons for doubting the validity of us- 
ing these potential forms in the analysis of the 
QGP-Hadron system. The final say in this mat- 
ter, however, rests with experimental data on 
the behaviour of C| in QGP. The overall picture 
presented by the various graphs, presented here, 
is the uniform variation of the velocity of sound 
over a range of temperatures including the ex- 
pected QGP-Hadron transition temperature of 
around 170 MeV, indicating a smooth roll-over 
of phases rather than a sharp transition. 
Note: 

To give a brief outhne of the origin of eq.(3), 
which wc have already given in detail [4], we re- 
call that the general Thomas-Fermi requirement 
on the number of particle states in configuration 
space is given by 

/ p{r)d^r = F{V{r)) (3a) 
Jo 

where p(r) is the coordinate-space density 
of states and F{V{r))is a functional of the 
coordinate-space potential V{r), both sides of 
(3a) being dimensionless. The momentum space 
counterpart of (3a) is 

pk 

/ p{k)dk^(f>{V{k)) (3b) 
Jo 

where in this definition of the density of states 
p{k), k is the energy variable and (f){V{k)) is 
the corresponding transform of F{V{r)). In this 
generalised form (3a), the potential may be ei- 
ther central as in atomic case, or a few body 
potential as in the nuclear and quark physics. 
In general to leading order the following forms 
are appropriate for (l){V{k)), namely 

3 

( Const . ) non—relat. 

(3c) 

[HV{k))]relat. ~ {C OUSt .) relat.[-V {k)f (3d) 

We can recover eq. (3) when we differentiate 
(3b) after substituting eq.(3d) in it with an ap- 



propriate choice of the constant dictated by the 
geometry of the system. 
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